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1. INTRODUCTION 



The Casimir effect [f] can be regarded as the change in the zero-point fluctuations due 
to nontrivial boundary conditions. Surveys of the effect are given, for instance, by Plunien 
et al. [2], Mostepanenko and Trunov [3], and Milton [4]. The recent "resource letter" of 
Lamoreaux [5] contains a wealth of references, although it is admittedly highly incomplete. 

In the past, the Casimir effect has been considered as a static effect. Growing interest 
in recent years has been drawn to the dynamical variant of the effect, meaning, in essence, 
that not only the geometrical configurations of the external boundaries (such as plates) but 
also their velocities play a physical role. Moore [6] is probably the first to have considered 
the dynamical Casimir effect. Examples of more recent references are [7] and [8]. 

The recent paper of Nagatani and Shigetomi [9] is an interesting development in this 
direction. These authors focused attention on the fact that if moving boundaries (mirrors) 
create radiation, the mirrors have to experience a reaction force. They proposed an effective 
theory for the back reaction of the dynamical Casimir effect in (1+1) dimensions for a scalar 
field, this theory being constructed by the background field method in the path integral 
formalism. In fact, they considered a kind of 2d quantum cosmology for describing the 
dynamical Casimir effect. 

In the present paper we show how to apply the effective action formalism, using the 
background field method, to formulate the dynamical Casimir effect in four dimensions in a 
convenient and elegant form. We are able to consider an arbitrary matter content (typically 
a grand unified theory or GUT) and present the dynamical Casimir effect as a kind of 
quantum cosmological model. Using the background field method, we treat the geometrical 
configuration of the boundaries classically, but consider the GUT in the interior region as a 
quantum object. 

In the next section we consider a GUT in a three-dimensional space, where the size 
of the space a(t) is a dynamical variable. Similarly to [9] we make use of the adiabatic 
approximation. Exploiting the conformal invariance of the theory we calculate the anomaly- 
induced effective action W. In the simplest case (a torus), W is given by Eq. (2.4). We 
consider the static Casimir energy in section 3, and show that, for the usual boundary 
conditions on the torus, the Casimir energy is attractive. In section 4 we start with the 
effective action, Eq. (4.1), for the dynamical case. Introducing a mass m associated with the 
scale factor a, with a corresponding kinetic energy in the low-velocity approximation equal 
to |md 2 (a phenomenological term), we then consider two cases. If the Casimir energy 
is attractive, we derive in Eq. (4.15) the time variation of a(t) for large values of t, the 
last term in the expression being a (special case of the) dynamical correction to the pure 
quasistatic Casimir result. For the perhaps less realistic repulsive case, the small and large 
time behavior of the Casimir behavior is extracted in Eqs. (4.24), (4.26). Numerical results 
in both cases are given in Section 5. The behavior of the scale factor in the two cases is 
shown in Figs. 1 and 3, while the dynamical stress on the torus is presented in Figs. 2 and 
4. In the Appendix we discuss the effects of alternative kinetic energy terms. 

2. THE EFFECTIVE ACTION 

Let us consider conformally invariant, massless matter in 4ci-dimensional space-time. The 
matter may correspond to some GUT (say, SU(5), SO (10), or any other alternative). We 
are interested first in the study of the static Casimir effect for such a theory when the field 
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is assumed to be bounded in a three-dimensional region. In other words, we are interested 
in a space having the form K 1 ® K 3 , where as K 3 one can take any manifold permitting an 
exact Casimir effect calculation. It can be S 3 , T 3 , S 1 <g> § 2 or any other compact manifold 
with a known spectrum of the d'Alembertian operator. We limit ourselves to T 3 or § 3 , for 
the sake of simplicity. 

Suppose now that our GUT lives in such a three-dimensional space, where the size of 
the space is a dynamical variable (moving mirror or moving universe). Hence, we will 
be interested in the dynamical Casimir effect in a three-dimensional region and the back- 
reaction from the induced radiation on the moving background geometry. We shall use the 
adiabatic approximation in this study. A great simplification comes from adopting a physical 
picture in which the Casimir effect is described as an effective action in curved spacetime 
(see [10] for an introduction). Here, spacetime is taken to be an expanding universe with 
topology M 1 ® § 3 or M 1 <g> T 3 . The corresponding metric is given by 

ds 2 = dt 2 - a 2 (t)ds 2 3 , (2.1) 

where ds 2 = dx 2 + dy 2 + dz 2 for T 3 (coordinates are restricted by all radii being equal), or 
the line element of a three-dimensional sphere S 3 . 

Let us calculate now the effective action for such a GUT. Using the fact that the theory 
is conformally invariant we may use the anomaly-induced effective action [11]: 

W = b J d^x^jFa 

+ b' J d 4 x^ {a[2a 2 + m»»Vfi v - □ + ^(V"£)(V„)]<7 + (G-^n R)a} 

~T2 X h b + y) J ^W^[£-66a-6(V^)(V^7)] 2 , (2.2) 

where our metric is presented in conformal form. Thus g^ v = e 2cr g^ u , a = In a(r)), i] is the 
conformal time, F is the square of the Weyl tensor, and G is the square of the Gauss-Bonnet 
invariant. Overbar quantities indicate that the calculation is made with g^ u . Further, 

" = iW (jV » + 6Ar " 2 + 12Ari) - 

" = "360^F (jY ° + lljYl/2 + 62iVl) - (2 ' 3) 

where iVo, Ni/ 2 , and Aq are the numbers of scalars, spinors, and vectors. For example, for 
M = 4 SU(N) super YM one gets [12] b = -b' = (N 2 - 1) [4(4tt) 2 ] -1 . We also adopt the 
scheme wherein the ^''-coefficient of the OR term in the conformal anomaly is zero. Being 
ambiguous, it does not influence the dynamics [12]. 

As the simplest case we consider henceforth a torus. Then 



W = Jd V [2b W"" - 2(6 + b')(a" + a' 2 ) 2 } . 



(2.4) 



This is a typical effective action for a GUT in a Friedman-Robertson- Walker (FRW) Universe 
of a special form. 
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3. THE STATIC CASIMIR ENERGY 



Let us briefly overview the static Casimir effect for a torus of side L (for more detail, see 
[13]). The Casimir energies associated with massless spin-j fields is 



(3.1) 



where we see the appearance of the characteristic minus sign associated with a closed Fermion 
loop. The frequency of each mode is given by 



<-=(x) 2 I>+«? > > 2 , » 



(ni, n 2 , n 3 ) 



Here g^ = 0, 1/2 depending on the field type chosen in R 1 ® T 3 . 

We use the p-dimensional Epstein zeta function Z p 
by the formula 



gi, g P 
hi, . . . , h. 



, ,u p 



(3.2) 



(s) defined for K s > 1 



gi, g P 
/ii, . . . , Zip 



(s) = J2'[(ni + gi) 2 + ... + (n p + g p ) 2 } 



21 ~P«/2 



x exp [27ri(ni/ii + ... + n p h p )] , 



(3.3) 



where gj and hi are real numbers, and the prime means omitting the term with (ni, n p ) = 
(— gi, — g p ) if all the gj are integers. For 5Rs < 1 the Epstein function is understood to 
be the analytic continuation of the right-hand side of Eq. (3.3). Defined in such a way, the 
Epstein zeta function obeys the functional equation 



gi, g P 
hi, . . . , h p 



(s) = n-^-^r Qp(l - *)) 



x exp [-27ri(gi^i + ... + g p h p )] Z p 



hi, . . . , h p 
~ gi) • • • i ~ gp 



(3.4) 



The function (3.3) is an entire function in the complex s plane except for the case when all 
hi are integers. In the latter case the function (3.3) has a simple pole at s = 1. 
Using Eq. (3.3) we have 



7T 



Z ^(-l) 2 ^3 



JJ) JJ) „U) 

61 g2 g3 





Taking into account the functional equation (3.4) one gets 



(3.5) 



JJ) JJ) JJ) 

61 62 63 




2tt 3 ' 











JJ) Jj) JJ) 
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(3.6) 



The Casimir energies (3.1) take the form 



'-TV, 



("1) 2J 

W^ 3 



-g? 



r U) 

52 



T (J) 
=3 



(3.7) 
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Finally the Casimir energy associated with a multiplet of fields characterized by the numbers 
No, N1/2, and N\ can be written as follows 



^1/2^3 



1 


N Z 3 








2tt 2 L 


(0) 

-gi " 


(0) 

-g 2 













(1/2) 
gl " 


-g^ 1/2) 


(1/2) 
-g 3 


G 



(0) 

-g 3 



-gl w 





g? } -g? 



(3.8) 



Thus, the static Casimir energy for a torus is proportional to c/L, where c is defined by the 
features of the GUT under consideration. Note that the sign of c is a priori unpredictable. 

If for all the fields of the theory we take the same boundary conditions (periodic or 
antiperiodic), the Z^s are all equal, and consequently the supersymmetry is not broken, 
and the Casimir energy is zero. On the other hand, if the different fields satisfy different 
types of boundary conditions, supersymmetry is broken and there is a static Casimir effect. 

For the latter situation, consider, as an illustration, the usual case of bosons satisfying pe- 
riodic boundary conditions and fermions satisfying antiperiodic boundary conditions on the 
torus. Then we require only two values. For the bosons, the Casimir energy is proportional 
to 








16.5323, 



(3-9) 



which value is given explicitly in Ref. [14]. For the fermions, the same reference gives the 
value 
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2 


2 





-3.86316. 



(3.10) 



So each term in Eq. (3.8) contributes a negative energy. Thus the net Casimir energy is 
attractive, 

£ = ~, c = 0.837537(iVo + iVi) + 0.195710iV 1/2 = 1.0332477V 1/2 , (3.11) 

since the number of fermions must be equal to the number of bosons. 1 

We should make the following general remarks concerning the physical interpretation of 
the calculation sketched here. Imposition of periodic boundary conditions at the boundaries 
of the field volume is a basic physical ingredient in expressions such as Eqs. (3.1), (3.2) for the 
Casimir energy. It is analogous to the imposition of perfect conducting boundary conditions, 
or more generally, electromagnetic boundary conditions, at the walls, when considering 
ordinary electrodynamics, for example within a spherical volume. The physical outcome of 



1 Any case with periodic boundary conditions in some directions and antiperiodic ones in others may be 
given in terms of the values given in Eqs. (3.9), (3.10) and the additional values 




\ 



0.689223, Z 3 





1 i 



-2.156887. 
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a calculation of this kind is the residual energy remaining when the influence of the local 
stresses is separated off. (Presumably, such stresses are absorbed in a kind of renormalization 
of physical parameters.) The field theoretical calculation is able to cope only with the cutoff 
independent part of the physical stress; the local cutoff dependent parts of the stress are 
automatically lost in the zeta-function regularization process. This is an important point 
whenever the result of the field theoretical calculation is to be compared with experiments. 

As a typical example of this sort, we may mention the calculation of the Casimir energy 
of a dilute dielectric ball. One may adopt a field theoretical viewpoint (cf., for instance, 
[15]), from which the Casimir energy is calculated as a cutoff independent, positive, ex- 
pression. More detailed considerations, using quantum mechanical perturbation theory [16] 
(cf. also [17]), or quantum statistical mechanics [18] show however how this expression is 
to be supplemented with attractive cutoff dependent parts. Such terms are presumably not 
observable. As for the cutoff independent term, agreement between the methods is found, 
so the situation is in this respect satisfactory. 



4. DYNAMICAL PROPERTIES 

Now let us turn to a simplified discussion of the dynamical Casimir effect. We here take 
into account that we have a dynamical radius a(t)L, a(t) being a dimensionless scale factor. 
Then, the total effective action is given as 

r = W-L J drja{rj)£, (4.1) 

where W is given by Eq. (2.4) and £ = —c/(aL) as displayed in Eq. (3.8). Because the 
action is dimensionless, the length L disappears from the calculation, and we have 

r = Jdrj [2b 1 aa"" - 2(6 + b')(a" + a' 2 ) 2 + c] . (4.2) 

This is a typical effective action to describe a quantum FRW Universe. 

In order to consider dynamical properties we add to the above effective action a phe- 
nomenological term, which has the form of a kinetic energy. We associate a mass m with 
the scale factor a, and take the corresponding kinetic energy to be given by |md 2 . Our es- 
sential idea is that the geometrical configuration of the space is treated classically and that 
the GUT field is a quantum object which induces the Casimir effect. One might in principle 
introduce other expressions for the kinetic energy, but this expression is clearly the simplest 
choice that one can make. The Newtonian form is moreover in correspondence with our use 
of the adiabatic approximation, meaning that \d(t)\ <C 1; cf. also the analogous argument 
in [9] in connection with the (1+1) dimensional case. 2 Introducing the physical time t via 
drj/dt = 1/a, we now write T as 



r = dt 



-ma 2 + 2b' In a ( a a 2 + 3 a da + d 2 a + da 2 ) 
2 



2{b + b') {h2 + a ' 6) \ C - 
a a 



(4.3) 



We consider other possibilities for the kinetic energy term in the Appendix. 
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where a = da/dt. 

From the variational equation ST /5a = we obtain, after some algebra, 3 

/ d 2 
ma -2b' [a 2 + 2— (ad) 
\ dt 2 

- 2 (6 + b') ( 2a a +4d a +3d 2 - 12— + 3^) + ^ = 0. (4.4) 
\ a a 1 J a 1 

It is remarkable that the the logarithm is absent in Eq. (4.4); there seems to be no reason 
a priori why this should be so. 

We limit ourselves to the M = 4 SU(N) super YM theory for which, as mentioned, 
{b + b') = 0. Then Eq. (4.4) simplifies to 

md + 26(2 ao + 4od + 3a 2 ) + = 0. (4.5) 

Both the terms involving b and c are dynamical, quantum mechanical, effects, which in 
dimensional terms are proportional to h. However, we will see that it is sensible (if c ^ 0) 
to regard the b term as a small correction to the Casimir-determined geometry. We denote 
the 6 = solution by ao(t); it satisfies the equation 

c 

md + — = 0, (4.6) 

On 



implying 



1 c 

-ma 2 , = h const. (4.7) 

2 do 



4.1. Attractive Casimir energy, c > 

If the Casimir energy is attractive, as actually realized in our illustrative calculation given 
in Sec. 3, see Eq. (3.11), we will assume, as boundary conditions, that a (t — > oo) = oo, 
cio(t — > oo) = 0. Then, the constant in Eq. (4.7) becomes equal to zero, and we get the 
Casimir solution 

2 / 9c \ ^ 

a (t) = At*, with A= [^) ■ ( 4 - 8 ) 

It is worth noticing here that the proportionality of ao(t) to t 2 ^ 3 is precisely the behaviour 
shown by the scale factor in the Einstein-de Sitter universe. This may be surprising at 
first sight, but does not seem to be so unreasonable after all, since the Einstein-de Sitter 
universe is flat, thus in correspondence with our neglect of Riemannian curvature terms in 
the formalism above. 



3 This and subsequent equations are dimensionally consistent if we restore dimensions: 

[a] = [t] = [m _1 ] = Length. 
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Now we turn to the solution a(t), taking into account the b correction. We shall limit 
ourselves to giving a perturbative solution, implying an expansion of a(t) around ao(t) 
assuming b to be small: 

a(t) = a (t) + ba 1 (t). (4.9) 

We consider only times for which the correction term is small: 

b ai /a < I. (4.10) 

Thus, we may expand the Casimir term in Eq. (4.5) as c/a 2 = (c/oq)(1 — 2ba\/ao). A first 
order expansion of the other terms in Eq. (4.5) then yields, when we take into account the 
Casimir solution (4.8), the inhomogeneous equation 

4 8A 2 _s 

The homogeneous version of Eq. (4.11) has solutions of the form t a , with a = 4/3 and 
a = —1/3. We write the independent solutions as 

f(t)=t~s, g(t)=t-s. (4.12) 

The Wronskian A between / and g is simple; A = fg — gf = —5/3. Writing for brevity 
the right hand side of Eq. (4.11) as r we then get, as the solution of the inhomogeneous 
equation, 



dt 



ai(t) = fit) (Ci-^J rgdt^j + g(t) (<0 2 + I J rf 

= ts d - — - r 2 + r 3 c 2 + — rs 

\ 15m / y 5m 

= Clt */z + c 2t -i/3 + i^! r 2 / 3 , (4.13) 

3m 

with Ci and C 2 being constants. As for the values of these constants, we have first to observe 
our restriction (4.10), which implies that 

b / 4A 2 \ 

J [c 1 ti+C 2 t- 1 + —t-i) «1. (4.14) 

If we require the perturbative approximation to be valid for large times, we must have 
Ci = 0. If we also set C 2 = 0, our perturbative solution becomes 

2 AA 2 2 

a(t) = At* + 6t~s, (4.15) 

3m 

which is only valid for large enough t, i.e., for 

bA 4 . 
— rs < 1. 4.16 

m 
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The static Casimir force is 
whereas the dynamical force is 

F dyn = md. (4.18) 

Substituting Eq. (4.15) into Eqs. (4.17) and (4.18), and observing the relation between c 
and A in Eq. (4.8), we get 

F dyn = F Cas (l - ^V 4/3 ) , (4.19) 

which shows that the dynamical force is the Casimir force modified by a small dynamical 
correction when the perturbative approximation is valid. 

Before we turn to a numerical solution of Eq. (4.5), we discuss the repulsive case. 



4.2. Repulsive Casimir energy, c < 



Now let us consider the case when c < 0, a repulsive Casimir energy. In this case we 
must take d |^oo 7^ 0. Let us write the 6 = equation (4.7) in the form 



2c 



c\ = 



m 



d = ±\/— + c 2 , 
a 



C 2 — Voo , Vqo — ao|t^oo- 



From (4.20) we obtain 



1 

c 2 



<7 () / — + c 2 In (2 v /c 2 a + 2Vc 2 ao + ci) 



±i + c 3 , 



where C3 is a further integration constant. We see from Eq. (4.20) that 

-2c 



mvi 



For long times, the solution behaves as 

ao(t) ~ y/c^t, t » 1. 
For short times, suppose ao approaches the minimum value (4.23); then 



c 3 = 



3/2 



ln2v^, 



and 



a (t) ~ -— - ^-t 2 , t < 1. 
c 2 4ci 



(4.20) 
(4.21) 

(4.22) 

(4.23) 

(4.24) 
(4.25) 

(4.26) 
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5. NUMERICAL SOLUTION AND DISCUSSION 



5.1. c>0 

Let us consider numerical solutions of dynamical equations (4.5) for SU(2) super Yang- 
Mills theory, for the attractive case. We suppose that the initial behavior of a(t) is given 
by the perturbative form Eq. (4.15). We may always set m — 1 since that amounts to using 
dimensionless variables for a and t. Let us take as an illustration 

AT = 2, c=l (A = 1.65096), t = 0.5, (5.1) 

For later times we integrate the exact equations numerically, starting with the initial condi- 
tions at t : 

a = 1.06744, a = 1.35019, a = -0.802706, a = 1.81582. (5.2) 

For those conditions we have a numerical solution for a(t) as shown in Fig. 1. For comparison 
we also show in the figure the unperturbed solution (4.8) due to the static Casimir force. 
It will be noticed that for large t there are significant deviations from the unperturbed 
solution, which must be due to C\ ^ in the perturbative solution (4.13). In fact, for the 
entire range of t — 0.1-80, the exact solution shown in Fig. 1 is roughly reproduced by 
Eq. (4.13) with C\ — C2 — —5. For the Casimir force we have the behavior as shown in 
Fig. 2. The exact solution has oscillations, but overall is close to the unperturbed solution. 
Not surprisingly, the Casimir energy dominates the force. Note also that for another choice 
of initial conditions one will find somewhat different behaviour. The essential property of 
the approximation under discussion is that there are always dynamical oscillations around 
the static Casimir force. 

5.2. c<0 

Next we consider numerical solutions of dynamical equations (4.5) for SU (2) super Yang- 
Mills theory when c < 0. We suppose that the initial behavior of a(t) is given by form 
Eq. (4.22). Let us take the illustrative values 

N = 2, ci = -l, c 2 = l, a \ t=0 = l, t = 0.1, (5.3) 

For later times we integrate the exact equations numerically, starting from the initial con- 
ditions at t : 

ao = 1.00241, d = 0.0490327, a = 0.497511, a = -0.049793. (5.4) 

Note that the perturbative values of these parameters, given from Eq. (4.26) are close to 
these: 

a = 1.0025, d = 0.05, d = 0.5, d = 0. (5.5) 

For those conditions we have a numerical solution for a(t) as shown in Fig. 3. For the 
Casimir force we have the behavior as shown in Fig. 4. For both cases the exact solution 
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has very small oscillations, but overall is close to the Casimir solution, and is accurately 
described by the limits of that solution, Eqs. (4.26) and (4.24). 

Thus, we have presented a formalism to describe the dynamical Casimir effect in the 
adiabatic approximation. It may be applied to an arbitrary GUT. Without any technical 
problems one can generalize the present consideration to any specific four-dimensional back- 
ground (we limited ourselves to a discussion of a toroidal FRW universe as providing the 
moving boundary conditions). But the limitations of our approach must be stressed: It 
would be extremely interesting to suggest new formulations of the dynamical Casimir effect 
beyond the adiabatic approximation. 
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APPENDIX A: ALTERNATIVE KINETIC ENERGY TERMS 

In the text, we introduced an ad hoc kinetic energy term into the action, referring to the 
change of scale with physical time, 



This is rather natural in the adiabatic context, where |d| <C 1, for then simple scaling 
properties obtain, as evidenced by the dimensional consistency of the resulting equations 
of motion when m has dimensions of mass. But we can not offer very strong arguments in 
its favor, in the absence of dynamical information. So, in this appendix we consider two 
alternatives, which provide somewhat different models for the dynamical evolution of the 
world. 

In the first, we suppose that the same kinetic energy should be integrated over conformal 
time, 4 



so that the Casimir evolution equation is, in place of Eq. (4.7), 



where we have dropped the subscript for simplicity, and written the constant of integration 
as k. The solution of this equation is very simple, 




(Al) 




(A2) 




a = o(0). 



(A4) 



4 In this case the parameter m is dimensionless. 
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If c > 0, we can set k = and obtain instead of the behavior exhibited in Eq. (4.8), a linear 
growth of the scale, 

a = a + \ —t. (A5) 
V m 

If c < 0, as before we cannot set the integration constant equal to zero; if we again choose 
the initial velocity to be zero, or a = —c/k, we get a result very like Eq. (4.26): 

a = -\ + ^ < A6) 

but now valid for all times. 

Perhaps a more natural possibility is to use the conformal time everywhere in the kinetic 
energy term, 5 

m ( da\ 2 m 



*?U =7 /<W ' (A7) 



The solution to the purely Casimir dynamical equation is 



^ad 2 = - + k, (A8) 



which is integrated to 



t = 7^\fj{l ^ {ka + c)3/2 ~ {ka ° + c)3/2] ~ 2c ^ {ka + 0)1/2 ~ {ka ° + c)V2] } • (A9) 

When c > again we can take k — > 0, which leads to the k — result 

a 2 = al + 2\—t. (A10) 
V m 

If c < and we choose again kao + c = 0, we obtain for short times 

a = -I + 2^ (2 ' (<<1 ' (A11) 

again very similar to Eq. (4.26). 

In Figs. 5 and 6 we show the effect of the inclusion of the dynamical b term for these 
kinetic energy structures. Qualitatively, the results do not depend much on whether the 
Casimir term is positive or negative. The example given for the first alternative kinetic 
energy is similar to the simple model result for c < shown in Fig. 3 except that the growth 
in t is quadratic rather than linear. The evolution for the second form of kinetic energy 
resembles the simple model result for c > shown in Fig. I. 
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FIG. 5: Dynamical behavior for a(t) for the first alternative kinetic energy term, Eq. (A2). Shown 
are the behaviors with m = 1, and initial condition a(0) = 1, evolving initially until t = 0.01 
according to Eq. (A4), with c = 1, k = 1 (solid line); and with c = — 1, k = 2 (dashed line). 
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FIG. 6: Dynamical behavior for a(t) for the second alternative kinetic energy term, Eq. (A7). 
Shown are the behaviors with m = 1, and initial condition a(0) = 1, evolving initially until 
t = 0.01 according to Eq. (A9), with c = 1, k = 1 (solid line); and with c = — 1, k = 2 (dashed 
line) . 
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